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The program objectives are fully defined in the original proposal entitled 


"A Research Program in Active Control/Aeroelasticity in the JIAFS at the NASA 
Langley Research Center" dated August 1, 1981. 

The research conducted by Dr. V. Mukhopadhyay during this report period and 
a summary of the overall research activities are described below: 


Development of Synthesis Methodology for Multifunctional 
Robust Aeroservoelastic System 

Introduction 

A synthesis methodology for multifunctional robust aeroservoelastic systems 
was developed. The development consisted of the following stages: 

1. Development of an universal diagram to determine phase and gain margins 
of a multi-input multi-output (MIMO) system using singular value based 
stability margin criteria. 

2. Determination of singular value gradients with respect to design para- 
meters and their application to improve stability margins of multiloop 
system. 

3. Application of constrained optimization techniques to synthesize a low 
order robust controller for a high order MIMO aeroservoelastic system 
while satisfying several design constraints on the dynamic loads and 
responses and stability margins at the plant input and output . 

System Description 

The MIMO feedback control system is described by the block diagram shown in 
Figure 1. The plant has several inputs and outputs denoted by vector U and Y' 
respectively. The root-mean square (RMS) dynamic loads and responses on which 
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upper-bound constraints may be applied are denoted by vector Y^. The maximum 
allowable value vector in Yqmax* The plant state-space matrices F, G u , H and an 
initial dynamic controller state-space matrices A, P, C are given. The initial 
controller may be obtained using Linear quadratic Gaussian (LGG) methodology, 
controller order reduction techniques or from a previous classical design. The 
initial closed-loop system must be stable. The expressions G(s) and K(s) repre- 
sent the plant and the controller transfer function matrices respectively, in the 
Laplace domain. The matrices R lT and Fy are so-called fictitious plant and 
measurement noise intensity matrices which can be used to improve the stability 
margins at the plant input (denoted by point (1)) and plant output (denoted by 
point (2)) for a full order LQG conroller design. 

In order to maintain certain phase and gain margins at the points ( 1 ) and/or 
(2) one must satisfy the conditions 

a ( I+KG ) > a D1 ((u) at input (1) 

a ( I+GK ) > a D2^ <4) ) at out P ut (2) 

where o denotes minimum singular value and is a measure of lower bound of the 

\ 

matrices in the parenthesis over the entire frequency range 0<a)<“. 

The inequalities (1) and (2) can be directly related to the desired phase 
and gain margins of a MIMO system using the universal diagram shown in Figure 2. 
Note that the singular value based stability margin criteria are sufficient con- 
ditions and may be conservative in nature. 

Problem Definition 

The controller synthesis scheme is posed as a constrained optimization prob- 
lem m an LQG framework and is shown m Figure 3. The cost function is a weight- 
ed sum of the expected value of steady state RMS plant outputs and controller 
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outputs. The constraints g 1 and 92 are scalars and are a cumulative measure of 
the difference between desired singular value shape and actual minimum singular 
value spectrum. The vector 93 represents all dynamic loads and design response 
constraints. The optimization algorithm uses the method of feasible directions 
to reduce all the constraints to zero or negative while minimizing the objective 
function J . 

Optimization Scheme 

The optimization scheme block diagram is shown m Figure 4. Starting from 
an initial stabilizing controller, the closed loop system matrices are computed. 

A set of Lyapunov equations are solved to obtain the objective function, the res- 
ponse constraints and their gradients with respect to the design variables. The 
singular values and their gradients are also obtained from analytical expres- 
sions. This information is used to update the controller using the method of 
feasible directions. The optimization algorithm CONMIN first attempts to satisfy 
all constraints while minimizing the objective function. 

Example Problem 

The synthesis method was applied to a two-input two-output system which 
represents a drone aircraft with a lateral attitude control system. A block 
diagram of the system is shown m Figure 5. The plant input position 1 is de- 
fined at the entry point to the elevon and rudder actuators denoted by and U 2 
in degrees units. The plant output position 2 is defined at the roll rate and 
yaw rate sensor outputs denoted by and Y 2 in degrees/second unit. The noise 
intensity matrices Ry and Ry can be used to improve the stability margins for the 
full order LQG controller. This is demonstrated in Figures 6 and 7 and Design 
Nos. 1 to 5 and are tabulated in Table I. In Design No. 1, for large value of 
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V the plant is made robust at plant output hut not at plant input as indicated 
by the singular value plots in Figure 6. By increasing R y as compared to Ry, the 
situation is gradually reversed as indicated by the Design Nos. 2, 3 and 4, in 
Figures 6 and 7. The singular value plot in Design No. 5 shows the beneficial 
effects of imposing more uncertainty in the unstable rudder channel. This design 
is used as a starting point for testing the present constrained optimization de- 
sign procedure, to irprove stability margins at both the plant input and output. 
In these designs the weighting matrices Q 1 = I and Qj = where I is appro- 

priate identity matrix. The contributions of U' and Y* to the objective function 
J are roughly equal. The noise intensity matrices P y and R v are set to zero. 

The system only contains a unit RMS white noise at the elevon actuator input. 

All the elements of matrices B and C are design variables and matrix A = (F+G U C- 
BH) as in a Kalman Filter. The minimum desired singular value g D = 0.45 is 
chosen as constraints on both cr(I+KG) and a(I+GK) at the plant input and outputs 
respectively. In computing the cumulative constraints, the frequency points are 
chosen uniformly in log scale as 50 divisions per decade. The Design No. 6, in 
Figure 7 shows the singular value spectrum after five iterations. The minimum 
a(I+KG) and ct(I+GK) are increased from 0.33 and 0.25 to 0.4 and 0.38 respec- 
tively. Thus the present constrained optimization procedure is able to improve 
the stability robustness at both input and output. However, the improvement is 
at the cost of some loss of high frequency attenuation at the plant input. The 
Design Nos. 7, 8 and 9 in Table I and Figure 8 show additional results for a 
third order controller using truncation, and optimization methods with the two 
types of constraints (A and B, as shown in Figure 3). 

Concise Statement of Research Accomplished 

A design methodology for multi-functional robust feedback controller was 
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developed. The method used constrained optimization techniques to improve stabi- 
lity margins of a multiloop system at both the plant input and output, while 
minimizing a standard L£XC performance index. The method can also be used to 
impose additional constraints on steady state and PMS responses. 
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Table - 1 Summary of Design Parameters and PMS Responses 
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* 0(A) - Optimization (Type A Constraint) 

0(B) - Optimization (Type B Constraint) 
T - Truncation of Design No. 6 





















FIG. 1 SYSTEM DESCRIPTION 
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FIG. 2 UNIVERSAL DIAGRAM FOR GAIN PHASE 

MARGIN EVALUATION 
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FIG OPTIMIZATION SCHEME BLOCK DIAGRAM 
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FIG. 5 EXAMPLE PROBLEM BLOCK DIAGRAM 
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FIG. 6 SINGULAR VALUE SHAPING 
BY NOISE ADJUSTMENT — FULL ORDER CONTROLLER 




FIG.? SINGULAR VALUE SHAPING BY NOISE ADJUSTMENT 
AND CONSTRAINED OPTIMIZATION - FULL ORDER 
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FIG. 8 SINGULAR VALUE SHAPING BY CONSTRAINED 
OPTIMIZATION — 3RD ORDER CONTROLLER 
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